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We consider a universe with a non-classical stringy topology that has fixed points. We concen- 
trate on the simplest example, an orbifold point, and study its observable imprints on the cosmic 
microwave background (CMB). We show that an orbifold preserves the Gaussian nature of the tem- 
perature fluctuations, yet modifies the angular correlation function. A direct signature of an orbifold 
is a single circle in the CMB that is invariant under rotation by 180°. Searching the 7- year ILC 
map of WMAP, we find one candidate circle and show that its high statistical significance is due to 
foreground contamination. We place a lower bound on the distance to an orbifold point at ~ 85% 
of the distance to the surface of last scattering. We show that due to galactic foregrounds, a more 
realistic bound is direction dependent and considerably lower. 



I. INTRODUCTION 

There is a strong experimental evidence that our universe is, to a good approximation, flat [1]. This, 
however, does not necessarily mean that the topology of the universe is M 3 . A non-trivial topology is a 
fascinating possibility that was investigated via its imprints on the cosmic microwave background (CMB) 
quite extensively (see e.g. [2-12]). So far the focus was on classical topologies [13]. Such topologies 
can be viewed as R 3 with some non-trivial identification that does not have a fixed point. For example, 
identifying x 3 with a; 3 + 27ri? we get R 2 x S 1 . This identification does not have a fixed point, and as a 
result R 2 x S 1 is flat. The extensive search failed to detect any sign of non-trivial classical topology. 

There is a good reason why the focus so far has been on classical topologies: If the identification has 
fixed points, then typically at the fixed points there is a curvature singularity and General Relativity 
breaks down. This is the sense in which these are non-classical topologies. This seems to suggest that we 
do not have tools to describe cosmology with non-classical topologies. However, a nice feature of string 
theory is that in many cases it resolves exactly these kinds of singularities. Roughly speaking, the way 
this comes about is that in string theory there are excitations that are confined to the fixed points and, 
together with the standard excitations that are free to propagate away from the fixed points, provide a 
consistent description of the physics everywhere, including at the fixed points (see e.g. [14]). From the 
point of view of string theory such topologies are as legitimate as classical topologies. 

This motivates us to initiate a study of the CMB imprints of such stringy topologies. Our focus here 
is on the simplest stringy topology - an orbifold point. The orbifold point is defined by the identification 

x-xo ~ -(x-Xo). (1) 

This identification has a single fixed point at xo, which is the location of the orbifold point. 

We do not expect to be able to detect cosmological imprints associated with the fixed point itself or 
the stringy excitations that are confined to it. However, if xo is within the observable universe then it is 
possible that (1) leaves a detectable imprint on the CMB sky. Part of the orbifold's imprint on the CMB 
is easy to illustrate. If the orbifold point is within the visible universe, there will be a special circle in 
the CMB that will be invariant under rotation by 180° (see Fig. 1). We refer to such a circle as a Self 
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FIG. 1: An illustration of the points on the last scattering surface (LSS) that are identified with other points on 
the LSS. As long as the orbifold is within the observable universe, these points form a thin circle. The distance to 
the orbifold ro (thick red line) and the opening angle of the circle a satisfy cos a = ro/r*. The phase separating 
each matching pair on the circle is n (two small black circles). 



Matching Circle (SMaC). In this work we study the signal to noise ratio (S/N) associated with a SMaC 
and more general imprints of an orbifold point on the CMB, and compare it to the data of the Internal 
Linear Combination (LLC) map of the Wilkinson Microwave Anisotropy Probe (WMAP). 



II. ORBIFOLD POINT 



In this section we discuss in some detail the orbifold point. In the first subsection we consider inflation 
in the presence of an orbifold point. We show that the orbifold does not induce non-Gaussianities. It does, 
however, modify the power spectrum in a specific way that breaks translation invariancc and isotropy. 
We use this power spectrum in the second subsection to calculate the CMB anisotropy two point function. 



A. Inflation with an Orbifold 



We assume the basic setup of slow-roll inflation in which the infiaton is a scalar field <fi. Then the 
orbifold identification (1) implies 

0(r ? ,x) = ^(r ? ,2x o -x), (2) 

where rj is conformal time. Quantum mechanically this means that 

[u( v , x), %(ri, x')] = i (<5(x - x') + <5(x + x' - 2x )) , (3) 

where u(r], x) = a(rf)(f){j), x) (here a(ry) is the scale factor) and 7r(r/, x) is the canonically conjugate variable. 
This yields the following commutation relation for the creation and annihilation operators 

= (2^) 3 (<5(k - k') + e - 2 ' lk - X0 <5(k + k')) , (4) 
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(and [cik) Ok'] = [ a k' a iL'] = 0)- The second term is due to the orbifold point and, as expected, it breaks 
both homogeneity and isotropy. 

Using these commutation relations we find that the theory is still Gaussian. Namely, only the two 
point functions do not vanish and yield 

(O|0(x)0(x')|O) = J ^ (e^C*"*') + e *-(*+*'-a*o)) P ^ k)j (5) 

where P^{k) is the standard inflaton power spectrum, H 2 /2k 3 , evaluated at horizon crossing. Hence for 
the metric perturbation, $, we get 

(* k *t) = (2^) 3 (<5(k - k') + e- 2ik ' x ° J(k + k')) P»(fc), (6) 

where P$(fc) is the standard power spectrum. 

We made two hidden assumptions in deriving (6). First, we assumed that the string scale is much 
smaller than the Hubble scale during inflation. Otherwise, the stringy modes that are confined to the 
orbifold point will affect the inflaton power spectrum that is fixed at the Hubble scale. Second, we 
assumed that the compactifaction scale is smaller than the Hubble distance during inflation. The reason 
we are forced to make this assumption is that the orbifold point breaks SUSY and as such, leads to 
instabilities. To avoid this issue we can consider an orbifold that acts also on the compactificd directions 
(that must exist in string theory) and that on large scales acts like (1). The simplest example is an 
orbifold point on t 3 x S 1 . If the radius of the S 1 is larger than the Hubble scale during inflation then (6) 
is a good approximation. These assumptions are often made and do not constrain the generality of (6) 
too much. 



atm = ' d fc$ kA,(fc)r/ m (k). (7) 



B. Angular Correlation Matrix 

Equipped with eq. (6) the calculation of the harmonic coefficients of the temperature anisotropy is 
straightforward (though cumbersome) since the a^ m 's are determined by $k in the following way 

J ^ 

Here A^(fc) = (k, 77*) is the (scalar) response function and 77* is the radius of the last scattering surface 
(LSS). In using the response function we are utilizing the "line-of-sight" approach [15], which takes into 
account the Sachs- Wolfe (SW) effect, the integrated SW effect on large scales and the Boltzmann physics 
of the coupled photon-baryon fluid on small scales. 

We use (6) to calculate the covariance matrix of the harmonic coefficients 

CimX'm' ( x o) = (o-lmO-e'm') ~ ?>U'$mm'Ce > + ^-Cimi'm' ( x o); (8) 



where 



Cf) = - / k 2 dkP^(k)\At(k)\ 2 (9) 



7T 

is the standard diagonal ACDM angular power spectrum and 

Ac tafmI (x„) = /+^ J d 3 fcp <& (fc)A,(fc)Ai,(fc) e - 24k - xo r^(k)^ m ,(k) (10) 

is the departure of the correlation function from ACDM due to an orbifold located at xo = rod. In order 
to evaluate (10) we use the familiar expansion of an exponent in spherical harmonics 



e ikr 



i'y/MM + m(kr)Y eo (k ■ r), (11) 
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where je(x) is the spherical Besscl function. We fix the coordinate system so that n points in the z 
direction, and get 

AC« w (ro,z) = ^ J k 2 dkP*{k)A e (k)A* e ,(k)A eme , m ,{2kr ), (12) 
where the matrix A£ m ^/ m / is given by 

oo 

A lm trn> (X) = ]T VMM" + Tjj t „ (x)G^jr'°. (13) 

e"=o 

The Gaunt integral Q is defined by 



d 2 kr £imi (k)y, 2 „ i2 (k)y, 3rn3 (k) 

(2£i + l)(2^ 2 + l)(2£ 3 + l) (l x It. h\ ( h h h 



Am \ I luii TO2 ms 



(14) 



and calculated using Wigncr's 3-j symbols. It is non-zero only if £1 + £ 2 + ^3 is even, mi + m 2 + m 3 = 
and \li — £j\ < < £i + £j. Therefore, the matrix At m £t m i is given by 

i+e' 

A imt >m'(x) = 5 mm , J2 (-l) m+{t+t '- n/2 VMM" + m"(x)Gu™°- (15) 

e"=\t-t'\ 

The resulting correlation matrix is real, proportional to 8 mm ' and symmetric with respect to both m O 
— m and I -H- £'. 

As a consistency check of this equation we take ro = 0. Then statistical isotropy is restored and the 
presence of the orbifold point implies that Ce — for all odd ts. Indeed, since ji>>(0) = <5f»o, we get 



A lM {Q) = S u ,S mm ,(-lY +m V^Gu n im0 

= SwS mm '(-lY, (16) 

which immediately gives the diagonal correlation matrix 

CWm'fo = 0) = SwSmm' [l + (-1)'] Cf\ (17) 

as expected. 



III. DETECTION OF AN ORBIFOLD IN THE CMB 



We have seen that the CMB temperature field in a universe with an orbifold point topology differs from 
that of a trivial topology in its two-point correlation functions. Our goal in this section is to calculate 
the S/N for the detection of the orbifold as a function of its distance ro. In the first subsection we use 
the full information in the CMB temperature correlation. However, as we explain below, this approach 
is not feasible. In the second subsection we study a more economical approach that searches for a SMaC 
in the CMB. 
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FIG. 2: The ideal S/N for orbifold detection, to leading order, (a) The S/N per multipole I for an orbifold 
located at ro = 0.5 (solid), 0.7 (dashed) and 1.1 ( dot- dashed) , in units of the radius of last scattering, (b) The 
accumulated S/N up to ^ max = 200 as a function of the distance ro. 



A. Ideal Signal-to-Noise for Detection 



As we have seen, the sole effect of the orbifold point is the deformation of the covariance matrix of 
the harmonic coefficients without inducing non-Gaussianities. Schematically, we got C or b = C^ ' + AC, 
where C^ ) is the standard ACDM covariance matrix. Since is diagonal, the S/N for detecting the 
orbifold is [16-18] 



Trfc^Cr 1 - l) +log (detC/dctC (0) ) . (18) 



N 

Much of the signal is at high-£ where this calculation is extremely computationally intensive. Hence we 
use the expression expanded to leading order in AC, 

S\ 2 v^/S 



if) -? #).■ < 19 > 



where 



The resulting expression depends only on the distance r to the orbifold point. 

In Fig. 2(a) we plot the S/N for each multipole t and for a few values of ro. We see that as expected 
the S/N is a decreasing function of ?' . This is also apparent from Fig. 2(b), where the accumulated S/N 
up to £ max = 200 is shown. 1 We can see that the S/N drops fast as the orbifold location approaches the 
last scattering surface. Once outside, it can no longer be detected. On the other hand, as long as the 



1 This scale corresponds to the 1° resolution of the ILC map. 
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orbifold is within the observable universe, S/N > 5. There are two contributions to the ideal S/N, one on 
large scales (£ < 10) and the other on intermediate scales (£ ~ 200). The dominant intermediate scales 
contribution is due to the SMaC and its typical scale is fixed by the thickness of the LSS. 

Ideally, we would use the correlation matrix to search for an orbifold signature in the CMB data. In 
comparing the observed correlations with those of the model, we would utilize as much of the available 
information as possible. However, we now see that this approach is not feasible, as most of the signal lies 
in small scales, requiring the analysis of large matrices containing high-^ data. Since the orbifold model 
is not isotropic, such an analysis includes the generation of a large correlation matrix for each orbifold 
distance rg as well as its rotation to all orientations with respect to the data. We therefore have to search 
for an orbifold signature in a way that will not use all the available signal, but will be computationally 
feasible. This is done in the next subsection by focusing on the SMaC only. 



B. A Self Matching Circle 

As explained in the introduction, if the orbifold point is within the visible universe, we expect to find 
a SMaC in the CMB. The opening angle a of the circle satisfies 

cos a = r /r*, (21) 

where r* is the radius of the LSS, as illustrated in Fig. 1. Every point on this circle is identified with 
another point on the circle, with a phase of n separating the two. Therefore, the CMB temperatures of 
the two points match. If the CMB temperature fluctuations were comprised of the SW signal alone, we 
would see an exact image of the surface of last scattering and the match would be perfect. There are, 
however, other effects that add noise to this match. Experience from the celebrated search for matching 
circles [6, 7, 9, 11, 12], relevant for classical topologies, suggests that overall the S/N should be quite 
high if the circle is large enough. That is, the S/N for detecting the orbifold is expected to decreases 
with rg. A larger ro means a smaller circle, and a smaller circle is harder to detect with high statistical 
significance. 

The score for detecting a SMaC is the following. Compare each circle in the data with itself, rotated 
by a phase of n. We therefore consider, for each pixel p and for each opening angle a, the score 

5 i \ ( T p{ a i4')T p {a,4> + 'K)) , . 

where T p (a, <fi) is the temperature at phase <f> along the circle centered at p and the angle brackets denote 
integration over <fi. We follow [7, 9, 11, 12] and Fourier transform the temperature profile along the circle 
as T p (a,<f)) = '^2 m T p ^ m (a) exp(im4>). We give appropriate weights to the different angular scales [9] and 
get 



= E m (-l) m ^l^.,m(« 



S p( a ) = ^ 7Z, 7~T72 ■ ( 23 ) 

This score equals unity for a perfect match and for a random circle its expectation value vanishes. 

In order to test if the CMB data contain a SMaC, we should examine all pixels p and all angular radii a. 
We use the HEALPix scheme [19] to search all directions using a grid with Aside = 256, corresponding to a 
resolution of ~ 0.25°. Following [7], when collecting data points along each circle, we linearly interpolate 
values at 2A s jd c points. For each a, we record S max {a), the highest value of S p (a) maximized over all 
pixels p. We search all angles in the range 20° < a < 90° with a resolution of 0.25°. 

We wish to verify that our score works as it should. To do that, we simulate CMB temperature 
fluctuation maps for a universe with an orbifold using the covariance matrix, calculated in §11 B. First, 
we get the response function A^?(fc,7?*) from CMBFAST [15] using CMBEASY [20], a computation 
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FIG. 3: The score S ma x(a) as a function of opening angle, for a simulated CMB sky with an orbifold located at 
ro — 0.5r*. A single spike appears at the correct a. The dashed line corresponds to the false detection level, 
discussed below. 

independent of the location of the orbifold. Then, we compute the covariance matrix C(ro) of an orbifold 
located in the z direction. We perform a Cholesky decomposition of the positive-definite matrix to get 
the triangular matrix L satisfying C = LL^ . We then draw a vector xg m of uncorrelated complex random 
Gaussian numbers of zero mean and unit variance, and compute the atm coefficients as a = Lx. For an 
orbifold in a different direction, the coefficients can be rotated using a Wigner rotation matrix. Finally, 
the temperature fluctuation is given by T(n) = J2e m a imYi m (iL) ■ We use £ max = 200 in our simulations. 
We plot S max (a) for a simulation with ro = 0.5r* in Fig. 3. We find that, as expected, there is a single 
spike, at the correct opening angle. 

C. Orbifold Detection Threshold 

We can use simulated sky maps to determine the maximal distance to an orbifold that could still be 
detected with our SMaC score, namely the detection threshold. An orbifold located farther away cannot 
be detected with our score as its signal is obscured by the statistical noise of ACDM . 

To estimate the noise, we simulate ACDM sky maps and calculate S max (a) for each of them and for 
each a. We then define a false detection level (FDL) corresponding to a level of 3cr above the median 
of the score of the random ACDM sky maps. 2 This is the dashed line shown in Figs. 3, 4 and 5(a). It 
is important to note that the 3<7 FDL we take has the meaning of local significance. By local, we mean 
the significance of a peak in S mlix (a) with respect to that specific opening angle a, and all directions. Of 
course, once a local peak is found, one needs to take into account the so called "look elsewhere" effect 
to estimate its true significance with respect to all radii and all directions. However, a peak with local 
significance lower than the FDL should not even be considered as a candidate for further analysis. 

We plot S mllx (a) for a few values of ro in Fig. 4(a). We can see that an orbifold which is located nearer 
to us (a SMaC with a larger radius) produces a higher spike in S max (a), and is therefore easier to detect. 
In addition, an orbifold located at a distance of ro = 0.85?% does not produce a peak in S ma x(a) that is 
significant enough to cross the FDL, setting the detection threshold to ~ 0.85?v 



2 Since S ma x(c0 is not normally distributed we cannot use the mean and standard deviation of the distribution. We 
therefore use the median and the percentile corresponding to 3cr instead. 
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FIG. 4: (a) The score 5 , ma x(a) for simulated orbifold sky maps at various distances ro. The dashed line is the 
FDL. As the orbifold is located farther, its peak is less significant, until at ro ~ 0.85r* it is obscured by the 
noise, (b) The score for a simulation with an orbifold at ro = 0.3r* and no contamination (solid) and for the 
same simulation with the data inside the KQ75 galactic mask replaced with a ACDM realization, once when the 
orbifold is in direction (Zi,6i) = (266°, -19°) (dotted) and once in (h,b 2 ) = (276°, -1°) (dot-dashed). 



The above analysis holds for a clean sky. A more realistic scenario, however, must take into account 
foreground contamination. Even when working with the ILC map, we should still be cautious about 
pixels that lie inside the galactic mask, as they can affect the matching circle and change the score. For 
instance, if an orbifold resided in the direction of the galactic north pole at a distance corresponding to, 
say, a ~ 80° , its SMaC would lie entirely inside the mask, risking undetectability even though according 
to Fig. 4(a) an orbifold at such a location should be easily detected. The mask can also affect smaller 
circles, with an arbitrary portion of the circle inside the mask. 

Since the mask is anisotropic, the FDL is also direction dependent. This in turn means that the 
detection threshold is direction dependent as well. We demonstrate this by considering two directions 
that have recently been found to be related to CMB anomalies and examining the effects of galactic 
foregrounds on an orbifold in each of these directions. The first is (h, b\) = (266°, —19°), normal to the 
reflection plane exhibiting large scale odd parity [21, 22], and the second is (l 2 , 62) = (276°, —1°), around 
which [23] found giant concentric rings. In order to simulate a map with foreground contamination, 
wc rotate an orbifold simulation so that the orbifold point is in the relevant direction, and replace the 
data inside the KQ75 galactic mask with those of a random ACDM realization. After contaminating 
the same orbifold simulation in orientations suitable for directions rii and 112, we search for the SMaC 
with our score. The results, for r = O.Sr* as an example, are plotted in Fig. 4(b), where it is clear that 
the orientation affects the level of contamination. Varying ro and repeating this analysis, we find the 
detection threshold for these directions. We get 0.35r* and 0.4r* for ni and 112, respectively, a significant 
change from the clean-sky isotropic 0.85^. 

Unfortunately, calculating an orientation-dependent FDL is not feasible. Therefore, we are left with 
our isotropic FDL and the naive estimate for the detection threshold. If the data of WMAP show any 
peaks in S max (a) that cross the FDL, we should check each candidate and see whether its orientation 
with respect to the mask is such that it could be heavily contaminated. On the other hand, there is also 
a chance that a SMaC does exist in the CMB and it is too contaminated for us to detect. 
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FIG. 5: (a) The score Smax(a) calculated on the 7- year ILC map. The dashed line is the FDL. The peak we 
analyze further is marked with a circle, (b) The match profile T(<p)T(cf> + ir) as a function of phase <j> for the 
circle centered at (I, b) = (275°,— 49°) with a = 85.75° of the ILC map [solid, left axis). The (smoothed) KQ75 
galactic mask is superimposed [dot-dashed, right axis). The entire match of the circle is due to a single pair of 
pixels, located inside the mask. 



IV. RESULTS 



The score <S' max (a) of the WMAP 7- year ILC map is shown in Fig. 5(a). It can be seen that one 
distinct peak crosses the FDL, at a = 85.75°. It is generated by a circle centered in the direction 
(/, b) = (275°, —49°). We examine it closer and check if some specific patches in that circle are responsible 
for its high score. For this sake, we revert back to our score in pixel space (22) and plot in Fig. 5(b) 
the match profile T(<fi)T(<fi + ir) as a function of the phase (j> along the circle. We find that the entire 
contribution to the score comes from a single pair of matching pixels. In addition, we also plot in Fig. 5(b) 
the profile of the KQ75 galactic mask, after smoothing it to 3°, along the same circle. The anomalous 
matching pair is well inside the galactic plane, where the mask is zero as indicated in Fig. 5(b). It includes 
the pixel at (I, b) = (0, 0), the center of the galactic bulk and the most heavily contaminated area in the 
map. We therefore conclude that this peak is caused by foreground contamination and does not indicate 
an orbifold. 

Since no other peak crosses the FDL, we conclude that the ILC data does not show a SMaC, with 
99.7% CL. 



V. DISCUSSION 



Testing the ILC map for a SMaC with our score, we have found no evidence that the universe has 
an orbifold topology. With the resolution of the ILC map, our score is not sensitive enough to find an 
orbifold that is located farther than ~ 0.857V This, however, does not mean that we can exclude the 
possibility that there is an orbifold point at tq < 0.85r„. We have demonstrated how foregrounds can 
affect the circles, changing their scores, and how important it is to mask the galactic plane. We showed 
that due to Galactic noise an orbifold point can even be, depending on direction, as close as tq = 0.35r* 
without being detected by our SMaC score. Data from Planck, soon to be available, will allow for a 
better analysis of self-matching circles. 

The possibility of finding support for string theory via CMB analysis is exciting. While the topology 
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of the universe remains a mystery, stringy topologies are viable candidates that should be thoroughly 
explored. We therefore intend to generalize our study here to other stringy topologies. 
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